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I. INTRODUCTION 

Finding out the reason of different behaviors of quan- 
tum magnets has attracted much interest in recent years. 
It is known that the interaction between atoms plays a 
major role. Recently, theoretical works are focused on 
the important role played by the orbital degree of free- 
dom in determining the magnetic properties of quantum 
magnets. The complex interplay among the spin and 
these types of degrees of freedom in quantum magnets, 
makes their phase diagram rich and induces various fas- 
cinating physical phenomena. Also, in some phenomena 
such as high-temperature superconductivity and colos- 
sal magnetic resistance, the orbital degree of freedom of 
electrons displays an important role. 

The experimental observations on Mott insulators are 
a realization of the effect of the orbital degree of freedom 
on the low-energy behavior of a system. A proposed tac- 
tic to investigate these systems is based on the so-called 
quantum compass model^. In fact, the quantum compass 
model is a very good candidate for explaining the low- 
temperature behavior of some Mott insulators. In this 
model, the orbital degrees of freedom are represented by 
(pseudo)spin-l/2 operators and coupled anisotropically 
in such a way to mimic the competition between orbital 
ordering in different directions. 

The two-dimensional quantum compass model is intro- 
duced as a realistic model to generate protected cubits^, 
and can play a role in the quantum information theory. 
This model is dual to studied model of superconduct- 
ing arrays^. It was shown that the eigenstates are two- 
fold degenerate and to be gapped^. The results from 
both spin wave study and exact diagonalization have 



suggested a first-order quantum phase transition in the 
ground state phase diagram^. The existing of the first- 
order phase transition is confirmed in two-dimensional 
quantum compass model^i^. Recently, it is found that 
in the ground state, most of the two-site spin correla- 
tions vanish and the two-dimer correlations exhibit the 
nontrivial hidden order—. Another examination of the 
2D quantum compass model showed that a second or- 
der phase transition can occur when the frustration of 
exchange interaction increases^. In addition, finite tem- 
perature properties have been studied. It is found that 
the low-temperature ordered phase with a thermal tran- 
sition corresponds to the 2D Ising universality class^iiiS 
and due to the dilution, the decrease of ordering temper- 
ature is much stronger than that in spin models^i. 

The ID version is known as the quantum compass 
chain (QCC) and represents one particular subclass of 
low-dimensional quantum magnets which poses interest- 
ing theoretical problemsi2r22. By mapping the model 
to a quantum Ising chain, an exact solution is obtained 
for the ground state energy and the complete excitation 
spectrunii^. It is shown that the QCC exhibits a first- 
order phase transition at Ji = between two disordered 
phases with opposite signs of certain local spin correla- 
tions. The model is also diagonalized exactly by a direct 
Jordan- Wigner transformationi^. The obtained results 
by latter approach, confirm the existence of the first- 
order phase transition in the ground state phase diagram. 
In a very interesting work, it is found that the reported 
first-order phase transition, in fact occurs at a multicrit- 
ical point where a line of the first-order transition meets 
with a line of the second-order transitions^. Based on 
a numerical analysisii, the first and second order quan- 
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turn phase transitions in the ground state phase diagram 
have been identified. By a detailed analysis of the nu- 
merical results on the spin structure factors, it is shown 
that the Neel and the stripe-antiferromagnet long-range 
orders exist in the ground state phase diagraroii. 

The effect of a transverse magnetic field on the QCC is 
studied by Ke-Wei Sun et al J^. They allowed the chang- 
ing a parameter that only provides moving on the first 
order critical line. Opening the energy gap of the sys- 
tem when a magnetic field applied is the result that they 
could achieve. In spite of the calculation of pseudo-spin 
correlation functions, fidelity susceptibility, the concur- 
rence, the block-block entanglement entropy, they didn't 
investigate the ground state phase diagram of the QCC 
in a transverse magnetic field. They only restricted their 
study on the first order critical line. In a very recent 
worki^, our group studied the presence QCC in an area of 
the ground state magnetic phase diagram where the odd 
couplings are antiferromagnetic and larger than even cou- 
plings. Using pseudo spin ladder operators indicated that 
the QCC in a transverse magnetic field can be mapped 
to the one dimensional Ising model in a transverse mag- 
netic field (ITF) . Analytical investigation of the effective 
ITF Hamiltonian predicts the occurrence of two quantum 
phase transitions by increasing magnetic field. 

In this paper we use the main idea of Perk^i to pro- 
vide analytical results on the effect of a transverse mag- 
netic field in whole range of the ground state phase dia- 
gram. Depending on values of the couplings, one or two 
quantum phase transitions can happen in finite critical 
fields. To make a clear picture of different long-range or- 
der phases, we do a very accurate numerical experiment. 
In particular we use the numerical Lanczos method to 
diagonalize numerically finite chains up to A'^ = 24 spins. 
Based on the exact diagonalization results we show that 
the spin-flop long-range order induces between two crit- 
ical magnetic fields. The spin-flop phase is denoted as a 
phase with the staggered magnetization perpendicular to 
the direction of the transverse magnetic fleld. 

The outline of the paper is as follows. In forthcoming 
section we introduce the results of the analytical fermion- 
ization studies. In the section III we discuss the model 
in the strong exchange couplings limit and derive the ef- 
fective spin chain Hamiltonian. In section IV, we present 
results of an accurate numerical experiment. Finally, we 
have discussed and summarized our results in section V. 



II. FERMIONIZATION 

The Hamiltonian of the QCC in a transverse magnetic 
field with N periodic sites is given by^^^ 



N/2 
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Here a^'^ are the Pauli operators on the zth site and 
Ji,J2,Li are the exchange couplings, h denotes the 
transverse magnetic field. In the absence of the trans- 
verse magnetic field, /i = 0, the first and second order 
critical lines denote with Ji/Li = and J2/L1 = 1, 
respectively^^. There are four gapped phases in the re- 
gions: (I.)Ji/Li > 0, J2/L1 > 1, (H.) Ji/Li > 0, 
J2/L1 < 1, (IH.) Ji/Li < 0, J2/L1 < 1, (IV.) Ji/Li < 0, 
J2/L1 > 1. There is a hidden order in the regions (I.) 
and (IV.). But two kind of magnetic long-range orders, 
a type of the Neel and the stripe-antiferromagnetic or- 
ders have been recognized in the regions (II.) and (III.) 
respectivelyii. 

A mapping of the spin-i operators onto fermi opera- 
tors by means of the Jordan- Wigner transformation was 
vastly used by many physicists who are interested to in- 
vestigate the Hamiltonian of the system by the aid of 
the spinless fermions. Because there are alternating links 
in the QCC we introduce two kinds of spinless fermion 
through the following Jordan- Wigner transitions^^ : 
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Using the above transformation, the QCC is mapped to 
a ID interacting spinless fermion system: 



nf - {^^^)Y^[aibi-a„b„] 
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Then by means of Fourier transformations, the mean field 
Hamiltonian is given by 



1 

hY,{A\A,+B\B^). 
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Finally, by diagonalizing the mean field Hamiltonian3i, 
the critical transverse fields obtain as 



he, = VJi{L + J2) 



he, = V-ML-J2)- (5) 

Surprizing is that, by applying the transverse magnetic 
field two quantum phase transitions can occur in the re- 
gion (I.). But in the regions (II.) and (III.), there is only 
one critical field. In the region (IV.), no quantum phase 
transition occurs by increasing the transverse magnetic 
field in a finite value of h. 



III. EFFECTIVE HAMILTONIAN 

In this section, we briefiy discuss the model ([T]) in lim- 
iting case of the exchanges ^ = ^ = ^ » 0. In this 
limit the model ([T]) can be mapped onto the effective spin 
chain Hamiltoniani^. 

In order to work easily, we rotated the system in the 
special sort that the field oriented in z direction. In the 
absence of the magnetic field, we consider one of the odd 
links for example the first link. For this dimer the quan- 
tum state is either the singlet state, \S), with the energy 
eigenvalue — J/2 or triplet state, |Ti), |T_i), |To) with 
energy eigenvalues of 0, 0, J/2 respectively. 

At /i = the ground state of a distinct dimer is \S). 
As the magnetic field h increases the energy of the triplet 
state |Ti) decreases and at h = J/2 forms together with 
the singlet state, a doublet of almost degenerate low en- 
ergy state. This degenerate state splits from the remain- 
ing high energy two triplet states. Thus, for a strong 
enough magnetic field we have a situation when the sin- 
glet \S) and triplet \Ti) states create a new effective spin 
T = 1/2 system. On the new singlet-triplet subspace 
and up to a constant, one can easily obtain the effective 
Hamiltoniani^ 
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The effective Hamiltonian describes the Ising chain in 
an effective transverse magnetic fLcldMh^^^ = h~ J . The 
critical effective fields for the Ising chain in a transverse 
magnetic field (ITF) are rigorously interesting, one of 
them is W^^^ = —L/2 and the other one is h'^JJ = L/2. 
Therefore critical fields of the main QCC can be obtained 
as 



he, = J -L/2 

he, = J + L/2. (7) 



It is clear that the hidden ordered gapped phase at 
h'^f^ < h'^JJ for the QCC corresponds to the nega- 
tively saturated magnetization phase for the effective 
ITF model. Whereas the intermediate gapped phase at 
h^JJ < h^fl < h^JJ for the QCC corresponds to the 
Neel phase for the effective ITF model. Finally, the re- 
gion h'^^^ > h'^J^f corresponds to the fully magnetization 
phase of the effective ITF model where the QCC is to- 
tally magnetized. In the next section, we present our 
numerical results obtained by a very accurate numerical 
experiment and confirm the mentioned quantum phase 
transitions. 



IV. NUMERICAL EXPERIMENT 

By doing an experiment, the validity of the suggested 
theoretical results on the effect of a transverse magnetic 
field on the QCC, can be determined. Since a real exper- 
iment cannot be done at zero temperature, the best way 
is doing a virtual numerical experiment. A very famous 
and accurate method in field of the numerical experi- 
ments is known as the Lanczos method. The Lanczos 
method and the related recursion methods2^i^, possibly 
with appropriate implementations, have emerged as one 
of the most important computational procedures, mainly 
when a few extreme eigenvalues are desired. However, the 
strong role of a numerical experiment to examine quan- 
tum phase transitions is not negligible. 

In this section, to explore the nature of the spec- 
trum and the quantum phase transition, we used Lanczos 
method to diagonalize numerically chains with length up 
to iV = 24 and different values of the exchanges. The 
energies of the few lowest eigenstates were obtained for 
chains with periodic boundary conditions. 

The first information that can be gained is the energy 
gap. In this way we have accessed to the energy gap 
which is recognized as the difference between the second 
exited state and the ground state energy in finite QCC. 
In Fig H] the energy gap is plotted versus transverse mag- 
netic field, h, for various length of chains. This figure 
contains four graphs of the energy gap in different sectors 
of the ground state phase diagram. The FigUJa) devoted 
to the energy gap of the system at the area (I.) of the 
phase diagram. In this figure the computed energy gap 
corresponding various length of chains [N = 12, 16, 20) 
and exchanges L = 1.0, Ji = 3.0 and J2 = 3.0 is plot- 
ted versus the transverse magnetic field. By increasing 
the size of the system the energy gap curve will approx- 
imate to the magnetic field axes. The magnetic fields in 
which the energy gap in the thermodynamic limit will be 
closed, are critical fields. By means of the phenomeno- 
logical renormalization group technique^^ we have found 
hci — 2.5 ± 0.1 and he2 = 3.5 ± 0.1 for critical fields in 
that the system undergos the quantum phase transition. 
In other words, in absence of the transverse magnetic 
field, the system is in a gapfull phase. By increasing the 
magnetic field, the gap of energy decreases and finally 



4 




FIG. 1: (Color online.) The energy gap as a function of the transverse magnetic filed h, for different chain lengths = f 2, f 6, 20 
and exchanges Li = 1.0, (a) Ji = J2 = 3.0, (b) Ji = 3.0, J2 = 0.5, (c) Ji = -3.0, J2 = 0.5 and (d) Ji = -3.0, J2 = 3.0 . 



will close at the first critical field hd ~ 2.5 ±0.1. Af- 
ter first critical field the gap of energy will be opened by 
growing magnetic field, while at the second critical field 
h — 3.5 ±0.1; again will be closed. For values of the mag- 
netic field further than h = 3.5 ±0.1 the gap of system 
will open and the system will go to a saturated phase. 
In conclusion, applying the transverse magnetic field to 
the gapped QCC in region (L), creates two new gap- 
full magnetic phases. The numerical critical fields with 
selected amount of exchanges L, Ji and J2 precisely cor- 
respond with the previous analytical critical fields (EqjS] 
and EqE]). 

In figures [IJb) and [IJc), the numerical results on 
the energy gap arc plotted for various length of chains. 



N = 12, 16, 20 and exchanges L = 1.0, Ji = 3.0, J2 = 0.5, 
in region (II.) and Ji = —3.0, J2 = 0.5, in region (III.). 
As it can be seen in FiglTJb) and (c), the qualitative effect 
of the transverse magnetic field in the regions (II.) and 
(III.) is equivalent. In the absence of the magnetic field, 
the system is gapped. Increasing the transverse mag- 
netic filed, the energy gap decreases and will be closed 
in he = 2.1 ± 0.1 [he = 1.2 ± 0.1) in regions (II.) and 
(III.) respectively. More increasing the magnetic field 
from he, the energy gap opens again and behaves almost 
linearly. The obtained numerical critical fields in these 
regions are in well agreement with the analytical results 
(Eql5]). Finally, Figlljd) allocated to the energy gap at 
the area (IV.) of the ground state phase diagram for chain 
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FIG. 2: (Color online.) The magnetization along the field , as a function of the transverse magnetic filed h, for different 
chain lengths N = 12, 16, 20, 24 and exchanges Li = 1.0, (a) Ji = J2 = 3.0, (b) Ji = 3.0, J2 = 0.5, (c) Ji = -3.0, J2 = 0.5 and 
(d) Ji = -3.0, J2 = 3.0. 



sizes N = 12, 16, 20, and exchanges L = 1.0, Ji = —3.0, 
J2 = 3.0. In this figure, is clearly seen that by increasing 
the transverse magnetic field, the energy spectrum of the 
QCC will remain gapfuU. There aren't any critical mag- 
netic fields. The system goes to a saturated magnetic 
phase for enough large values of the transverse magnetic 
field. 

So far the paper was determined that the transverse 
field, depending on the values of the exchanges can in- 
duces two or one gapped phases in the ground state 
phase diagram of the QCC. Subsequently, a very im- 
portant question that precedes is: "What are the mag- 
netic long-range ordered phases in these gapped regions 



of the QCC?" A very clear answer can provides, do- 
ing the numerical Lanczos experiment. To recognize the 
different magnetic phases induced by transverse field, 
h, we implemented the algorithm for finite-size chains 
{N = 12, 16, 20, 24) to calculate the order parameters 
and the various spin correlation functions. 

An insight into the nature of different phases can be 
obtained by studying the magnetization process. The 
magnetization along the transverse field axis is defined 
as 

1 ^ 

My^-J2(^s\ay\Gs), (8) 

i=l 
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where the notation {Gs\...\Gs) represents the ground 
state expectation value. In Figl^J we have plotted the 
transverse magnetization, M^, versus h, for different 
chain sizes N = 12, 16, 20, 24 and different values of the 
exchanges that covered all regions in the ground state 
phase diagram. We did not find any size effect on the 
numerical results of the transverse magnetization. As is 
seen in Figj^J^a), due to the profound effect of quantum 
fluctuations, the transverse magnetization remains small 
but finite for < /i < and reaches zero aX h = 0. This 
behavior shows that in the gapped phase with hidden or- 
der in region (I.), the transverse magnetization appears 
only at a finite critical value of the transverse field, hc^ . 
This phenomenon is also reported for a two-leg ladder— 
and an alternating spin chaiii^^. For h > hc^ , transverse 
magnetization increases almost linearly with increasing 
transverse field. At the second critical field hc2, due to 
the quantum fluctuations the magnetization do not satu- 
rate and reaching saturation asymptotically in the limit 
of infinite transverse field. 

Numerical results on the magnetization function pre- 
sented in figures El^b) and [2{c) are computed for vari- 
ous length of chains {N = 12, 16, 20, 24) and exchanges 
L = 1.0, Ji = 3.0, J2 = 0.5, in region (II.) and Ji = -3.0, 
J2 = 0.5, in region (III.). As the same as the energy gap 
function, the qualitative effect of the transverse magnetic 
field in the regions (II.) and (III.) is equivalent. In the 
absence of the transverse field, the transverse magneti- 
zation is zero. As soon as the transverse fleld is applied, 
the transverse magnetization start to increase. In the re- 
gion (III.), the transverse magnetization increases faster 
than the region(II.). It is also completely clear that the 
system at the critical fleld he do not saturate. The sat- 
uration will be happened for very larger values of the 
transverse field. Figl^fd) shows the numerical results in 
the region (IV.) of the ground state phase diagram for 
chain sizes N = 12,16,20,24, and exchanges L = 1.0, 
Ji = —3.0, J2 = 3.0. As soon as the transverse field is 
applied, the transverse magnetization starts to increase 
from zero. The transverse magnetization grows smoothly 
by increasing the transverse field and due to profound ef- 
fect of the quantum fluctuations, it will be saturated in 
inflnite transverse fleld. By comparing flgures [2{a) and 
[SJd), one can flnd that the hidden order in the region 
(I.) must be completely different from the hidden order 
phase in the region (IV.). The hidden order phase in the 
region (I.), will not surrender versus the transverse fleld, 
but the hidden order phase in the region (IV.) as soon 
as the transverse fleld applied will surrender versus fleld 
and spins start to aline the fleld direction. 

To flnd the long-range magnetic order of the ground 
state of the system, we start our consideration with the 
spin-spin correlation function deflned by 



W-"(n) = (aX-hJ (« = 2;,y,z), (9) 



and the spin structure factor at momentum q deflned by 
^""(g) = W°"'{n)exp{iqn). (10) 

n=l 

It is known that the spin structure factor give us a deep 
insight into the characteristics of the ground stateii. 

By examining the energy gap we found that a new mag- 
netic phase is induced by applying the transverse fleld in 
the region (I.) between hc^ < h < hc^. To determine how 
to long-range order there is, in Figl^Ka) we have plotted 
Wf ^ as a function of n for different values of the trans- 
verse fleld. The selected values of the transverse fleld 
h = 0.7,3.0,4.2 in the flgure cover all phases in the re- 
gion (I.). It can be seen from this flgure (Figl^Ja)) that 
the z component of the spins on odd sites is pointed in 
the same direction with the erf and others (on even sites) 
are pointed in opposite direction for a value of h ~ 3.0. 
This is an indication for the Necl ordering in the interme- 
diate region hc^ < h < hc^. This is very interesting that 
the numerical results suggest that a different novel phase 
can be induced by increasing the transverse fleld which 
is known as the spin-flop phase. For any value of the 
transverse fleld the Lanczos results lead to the staggered 

magnetization, = Y^^=i{~^Y '^j = 0^ since the 
ground state is degenerate and in a finite system no sym- 
metry breaking happens. But, the /i-dependency of the 
spin structure factor, ^^^(Tr), is qualitatively the same 
as the staggered magnetization, M|(. However the spin- 
spin correlation function diverges in the ordered phase as 
N — > 00. We computed the spin structure factor S^^{Tr) 
and in Figl^b), wc have plotted S^^-{q = Tr)/{N ~ 1) as a 
function of h for the chain lengths N ~ 12, 16, 20, 24 and 
exchanges L = 1.0, Ji = 3.0, J2 = 3.0. It can be seen 
that in the intermediate region hci < h < hc^ , the ground 
state of the system is in the Neel phase. There is not the 
Necl ordering along z axis in regions h < hc^ = 2.5 ± 0.1 
and h > = 3.5 ±0.1. To check the existence of the 
Neel order in the thermodynamic limit N — > 00 of the 
system in the intermediate region hc^ < h < hc.^ , we have 
plotted in Fig[31[c) the TV dependence of S^^{Tr) for differ- 
ent values of transverse fleld. Increasing behavior of the 
spin structure factor in the sector hc^ < h < hc^ , shows 
that the spin-flop phase is a true long-range order. 

On the other hand the behavior of the energy gap 
versus the magnetic field in the regions (II.) and (III.), 
suggested a quantum phase transition at the critical 
field he- In a very recent workii, Neel and stripe- 
antiferromagnetic long-range orders, in the absence of 
the field, have been reported in the regions (II.) and (III.) 
respectively. The " stripe-antiferromagnetic phase" is de- 
fine as a phase with the opposite magnetization along 
the z axis on the odd bonds. The order parameter of the 
stripe-antiferromagnetic phase is defined as^i 

^sp = ^(E(-l)'K-l +^2,))- (11) 
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FIG. 3: (a) The spin-spin correlation Wf^ as a function of 
n, for chain length A'' = 20 and exchanges Li = 1.0, Ji — 
J2 ~ 3.0. (b) The spin structure factor S^^(n)/{N — 1) as a 
function of the transverse field h for different chain lengths 
N = 12,16,20,24 and exchanges Li = 1.0, Ji = J2 = 3.0. 
(c) The spin structure factor S^^{tt) plotted as a function of 
the chain length A'^, for different values of the field in the 
intermediate region. 



FIG. 4: (Color online. )The spin structure factor S'"'{-k)/{N- 
1) as a function of the transverse field h for different chain 
lengths iV = 12,16,20,24 and exchanges Li = 1.0, Ji = 
3.0, J2 = 0.5. 

We computed the spm structure factor S'^^(7r) and the 
correlation function of the stripc-antiferromagnetic order 
parameter given by 



Af/2-1 



2j+2n 



))• 

(12) 



In FigHl we have plotted S^^{tt)/{N — 1) as a func- 
tion of h for the chain lengths N = 12, 16, 20, 24 and 
exchanges L = 1.0, Ji = 3.0, J2 ~ 0.5. Clearly be seen 
that the Neel order remains in the presence of the field 
up to the critical field he = 2.1 ± 0.1. Overlapping of the 
numerical results in the region h < he = 2.0 ± 0.1 shows 
a divergent behavior of the function S^^ (tt) by increasing 
the size of chain N. This justifies that the Neel ordering 
along the z axis is true long-range order in the region 
h < he — 2.1±0.1 of the ground state phase diagram. 
Induced quantum fiuctuations by increasing h from zero, 
decreases the staggered magnetization from the almost 
saturation value. We checked our numerical results and 
found that there is not the Neel ordering along z axis in 
the region h > he- 

The final insight into the nature of different phases 
be related to the effect of the field in the region (III.) 
Which can be obtained by studying the correlation func- 
tion x^^. In Figinija), we have plotted x^^/(-^/2 — 1) 
as a function of the field h for different values of the 
chain length N = 12,16,20,24 and exchanges L = 1.0, 
Ji = 3.0, J2 = —0.5. It is completely clear, that 
there is stripe-antiferromagnetic ordering along the z 
axis in the region h < he = 1.2±0.1. The induced 
quantum fiuctuations by increasing h from zero, de- 
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FIG. 5: (Color online.) (a) The correlation function of the stripe antiferromagnetic order parameter y^^l(Nl2 — V) as a function 
of the transverse field h for different chain lengths A'' = 12, 16, 20, 24 and exchanges L\ — 1.0, Ji = —3.0, J2 — 0.5. (b) The 
correlation function versus the chain length N for different values of the field and exchanges Li = 1.0, Ji = —3.0, J2 = 0.5. 



creases the stripe-antiferromagnetic order from almost 
saturation value. Again, to check the existence of 
the stripe-antiferromagnetic order in the thermodynamic 
limit N — > 00 of the system in the region h < he = 
1.2 ± 0.1, we have plotted in FiglSjb) the N dependence 
of for different values of field h < he = 1.2 ±0.1. As 
is seen from this figure in the region h < he = 1.2±0.1 
there is a divergent behavior which shows that the stripe- 
antiferromagnetic order is true long-range order. By in- 
vestigating the dependence of for different val- 
ues of Ji/Li, we found that there is no long-range 
stripe-antiferromagnetic order for fields larger than the 
he = 1.2 ±0.1. 

V. CONCLUSION 

In this paper we considered the ID quantum compass 
model (QCC) and built a clear picture of different mag- 
netic phases induced by an external transverse magnetic 
field at zero temperature. This picture is made step 
by step. At the first step, using the analytical spinless 
fermion approach wc found that, depending on the val- 
ues of the couplings, by increasing the transverse mag- 
netic field, one or two quantum phase transitions occur 
in the ground state magnetic phase diagram of the QCC. 
These quantum phase transitions belong to the univer- 
sality class of the commensurate-incommensurate phase 
transition. 

In the second step, we did an accurate numerical ex- 
periment using the Lanczos method. This part of our 



work, helped us to increasing the brightness of the ground 
state phase diagram picture. We have implemented the 
Lanczos method to numerically diagonalize finite chains. 
Using the exact diagonalization results, first we have cal- 
culated the energy gap. In complete agreement with our 
analytical results, we showed that depending on the val- 
ues of the couplings, by increasing the transverse mag- 
netic field, the energy gap closed at one or two critical 
fields. Then we also studied the magnetization process in 
all regions of the phase diagram. By comparing the mag- 
netization curves, we showed that the suggested zero-field 
hidden ordered phases have different behaviors in pres- 
ence of a transverse field. The hidden order phase in 
the region (I.), will not surrender versus the TF, but the 
hidden order phase in the region (IV.) as soon as the TF 
applied will surrender versus field and spins start to aline 
in the field direction. Finally, to distinguish the kind of 
magnetic orders wc plotted the spin-spin correlation func- 
tions and found that the spin-flop long-range order can 
be induced in a sector of the ground state magnetic phase 
diagram. 
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